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Abstract 

Recent investigations of the Dyson-Scliwinger equations and Monte- 
Carlo lattice calculations resulted in a coherent description of the fully 
dressed gluon, ghost and quark propagators in Landau gauge QCD. In the 
Dyson-Schwinger approach the infrared behaviour of these propagators is 
determined analytically. For finite spacelike momenta the gluon, ghost 
and quark propagators are compared to available corresponding results 
of lattice Monte-Carlo calculations. For all three propagators an almost 
quantitative agreement is found. These results for the non-perturbative 
propagators allow an analytical verification of the Kugo-Ojima confine- 
ment criterion. Our numerical analysis clearly reveals positivity violation 
for the gluon propagator generated by a cut in the complex momentum 
plane. The non-perturbative strong running coupling resulting from these 
propagators possesses an infrared fixed point. The quark propagator ob- 
tained from quenched and unquenched calculations displays dynamical 
chiral symmetry breaking with quark masses close to "phenomenological" 
and lattice values. We confirm that linear extrapolations of the quark 
propagator for different bare masses to the chiral limit are inaccurate. 



*Based on an invited talk given by R.A. at the Second International Conference "Nuclear 
and Particle physics with CEBAF at Jefferson Laboratory, Dubrovnik, Croatia, 26 — 31 May, 
2003 
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1 Some Aspects of Confinement 



Data taken at the Jefferson Laboratory over the last few years have uncovered 
many unexpected properties of hadrons. It seems that CEBAF operates in a 
very interesting energy range, and it is foreseeable that after its upgrade many 
more highly interesting results will be obtained. More or less all of these in- 
vestigations aim at an understanding of the structure of hadrons in terms of 
the underlying degrees of freedom, the quarks and gluons. To bring our current 
theory of strong interactions. Quantum Chromo Dynamics (QCD), into agree- 
ment with observations the hypothesis of confinement is needed. Even thirty 
years after the formulation of QCD a detailed understanding of the confinement 
mechanism(s) is still lacking. 

Experimental data obtained at the Jefferson Laboratory pose a challenge to 
theoretical physicists working in this field: In principle the requirement is to 
understand hadronic properties at intermediate and large momentum transfers 
in terms of QCD degrees of freedom. An adequate theoretical approach for 
such investigations has to be deeply rooted in Quantum Field Theory to ac- 
commodate a description of confinement, it has to be Poincare-covariant to be 
applicable at the employed momentum transfers, and, last but not least, it has 
to be manageable for quite complicated hadronic form factors and reactions. 

Non-perturbative calculations of QCD correlation functions do fulfill all these 
requirements. Studies of their equations of motion, the Dyson-Schwinger equa- 
tions (DSEs), have the potential to provide a successful phenomenology of 
hadrons in terms of quarks and gluons, for recent reviews see e.g. |2 
Within this approach dynamical breaking of chiral symmetry is obtained from 
first principles 0]. The description of meson properties and reactions based 
on their Bethe-Salpeter amplitudes has been quite successful, see ref. PP and 
references therein, and it seems viable that covariant Faddeev amplitudes will 
allow for a computation of baryonic reactions as well, see e.g. On the other 
hand, the results of ref. |H] make clear that an understanding of confinement 
and its implications onto the analytic structures of QCD Green's functions is a 
necessary prerequisite for further progress in this direction. 

The phenomenon of confinement is truely non-perturbative in nature: There 
is a physical scale associated with it, and we know from asymptotic freedom of 
QCD that such scales are non-analytic in the coupling. Futhermore, as we will 
see shortly, QCD Green's functions exhibit infrared singularities related to con- 
finement. Lattice calculations are very helpful because they provide a rigorous 
non-perturbative method, nevertheless, there is definite need for continuum- 
based methods. In the following we will discuss how the infrared behaviour of 
QCD propagators can be determined analytically, and what the corresponding 
results tell us about confinement in the covariant gauge. 

Let us take a step back and discuss Quantum Electro Dynamics (QED) first. 
In QED in the covariant gauge the electromagnetic field can be decomposed 
into transverse, longitudinal and time-like photons, but the latter two are never 
observed. From a mathematical point of view this can be understood from 
the representations of the Poincare group: Massless particles have only two 
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possible polarizations. This apparent contradiction is resolved by the fact that 
time-like and longitudinal photons cancel exactly in the S'- matrix 0. We can 
interpret this also otherwise: The time-like photon, being unphysical due to the 
Minkowski metric from the very beginning, "confines" the longitudinal photon. 

In QCD cancelations of unphysical degrees of freedom in the ^'-matrix also 
occur but are more complicated due to the self-interaction of the gluons. One 
obtains e.g. amplitudes for the scattering of two transverse into one transverse 
and one longitudinal gluons to order Og. A consistent quantum formulation in 
a functional integral approach leads to the introduction of ghost fields 7 . To 
order a| a ghost loop then cancels all gluon loops which describe scattering of 
transverse to longitudinal gluons. The proof of this cancelation to all orders in 
perturbation theory has been possible by employing the BRS symmetry of the 
covariantly gauge fixed theory [5]. At this point one has achieved a consistent 
quantization. If one wants to describe confinement of coloured states in a similar 
efficacious way one has to go only one slight step further: One has to require 
that only BRS singlets are allowed as physical states pilTn|. 

The Kugo-Ojima confinement scenario [HI EI describes a mechanism by 
which the physical state space contains only colourless states. The coloured 
states are not BRS singlets and therefore do not appear in S'-matrix elements: 
They are confined. Transverse gluons are BRS-non-singlets states with gluon- 
ghost, gluon-antighost and gluon-ghost-antighost states in the same multiplet. 
Gluon confinement then occurs as destructive interference between these states. 
In Landau gauge a sufficient criterion for this type of confinement to occur is 
given by the infrared behaviour of the ghost propagator: If it is more singu- 
lar than a simple pole the Kugo-Ojima confinement criterion is fulfilled. It is 
important to note that the general properties of the ghost DSE and one ad- 
ditional assumption, namely that QCD Green's functions can be expanded in 
asymptotic series in the infrared, allow to prove this version of the Kugo-Ojima 
confinement criterion [111 I12| . 

2 Gluon Propagator and Running Coupling 

Recently strong arguments have been provided in favor of infrared dominance of 
the gauge fixing part of the QCD action in the covariant gauge ^13^ . The related 
infrared dominance of ghost loops also occurs in truncation schemes of DSEs 
being self-consistent at the level of two-point functions [T^ . These schemes have 
been refined and generalized jJS] and allowed then to solve the coupled set of 
DSEs for the ghost, gluon and quark propagators 0. 

In Landau gauge these momentum-space propagators Dg{p), D^i,(p) and 
S(j)) renormalized at a scale can be generically written as 
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(1) 




(2) 



3 



^^'^ ^ -^;}A(p^/l2) + B(p2^;,2) _,^+M(p2)- ^-^^ 

Two renormalisation scale independent combinations build from these functions 
will be important for the further discussion: M(p2) = B{p'^ , ^'^)/A{p'^ , fi^) de- 
notes the quark mass function, and a non-perturbative definition of the running 
coupling, as{p^) = as(^^) G'^{p'^,^^) Z{p'^,y?), is possible due to the non- 
renormalisation of the ghost-gluon vertex in Landau gauge |14j . 

Employing asymptotic expansions for the propagators at small momenta the 
ghost and gluon equations can be solved analytically. One finds simple power 
laws, 

z{p\^^') - (pVA^^)^^ G{p\^?) ^ {p^li,^Y\ (4) 

for the gluon and ghost dressing function with exponents related to each other. 
Hereby k is an irrational number, k = (93 - v'2101)/98 « 0.595 [T21[Tn|. The 
product G'^ij? ,\j?^Z{j? goes to a constant in the infrared. Correspondingly 
we find an infrared fixed point for the running coupling, 

4^ r(3-2K)r(3 + «)r(i + A.) ^ 
= r2(2-K)r(2.) " 2.972 

for the gauge group SU(3). This result depends slightly on the employed trun- 
cation scheme. Infrared dominance of the gauge fixing part of the QCD action 
|13| implies infrared dominance of ghosts which in turn can be used to show jl2| 
that as(0) depends only weakly on the dressing of the ghost-gluon vertex and 
not at all on other vertex functions. 

The running coupling as it results from numerical solutions for the gluon, 
ghost and quark propagators can be quite accurately fitted by the relatively 
simple function 0] 

, 2^ _ «5(0) I 4^ P'/Agep / 1 1 \ . . 

^ 1 + p2/A2,„ ^ /3o 1 + p2/A2,„ In(p2/A2,,) p2/A2^^ _ 1 ; ^ ^ 

with /?o — (lliVc — 2iV/)/3. Note that, following ref. the Landau pole has 
been subtracted. The scale Aqcd is hereby determined by fixing the running 
coupling at a certain scale, e.g. as{M'^) = 0.118. 



3 Quark Propagator 

In the quark DSE as well as in the quark loop of the gluon equation the quark- 
gluon vertex enters. It has proven successful ^ to assume that the quark- gluon 
vertex factorizes, 

(g, fc) = xc'^' q, k) w-'^'^' (p, q,k), (6) 

with p and q denoting the quark momenta and k the gluon momentum. The 
non-Abelian factor vt^^°'"=' multiplies an Abelian part VJ''"^', which carries the 
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Figure 1: The quark propagator functions in quenched approximation as well 
as for three massless flavours compared to the lattice data 

tensor structure of the vertex. For the latter we choose a construction TH' used 
widely in QED. 

The Slavnov-Taylor identity for the quark-gluon vertex implies that the 
non-abelian part W^°'''^''{p, q, k) has to contain factors of the ghost renormaliza- 
tion function G{k^). Due to the infrared singularity of the latter the effective 
low-energy quark-quark interaction is enhanced as compared to the interaction 
generated by the exchange of an infrared suppressed gluon. Therefore the effec- 
tive kernel of the quark DSE contains an integrable infrared singularity. Further 
constraints imposed on W^"''"^''{p, q, k) are such that the quark mass function is, 
as required from general principles, independent of the renormalization point 
and the one-loop anomalous dimensions of all propagators are reproduced. 

In Fig. n we compare our results for the quark propagator in quenched ap- 
proximation as well as for three massless flavours with lattice data JHl- These 
results nicely agree with the one from the lattice. Furthermore, for the consid- 
ered number of flavours the quenched approximation works well. 

4 Spectral Properties of the Gluon Propagator 

The infrared exponent k is an irrational number, and thus the gluon propagator 
possesses a cut on the negative real axis. It is possible to fit the solution for 
the gluon propagator quite accurately without introducing further singularities 
in the complex p^ plane. The fit to the gluon renormalization function |21| 



is shown in Fig. |21 Hereby w is a normalization parameter, and 7 = (— 13iVc -I- 
4:Nf)/{22Nc — 4:Nf) is the one-loop value for the anomalous dimension of the 
gluon propagator. The corresponding discontinuity along the cut vanishes for 
p^ 0^ , diverges to +00 at = — Aggp and goes to zero for p^ 00. 
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Figure 2: The gluon propagator (left diagram) and its Fourier-transform (right 
diagram) compared to the fit, Eq.(|Zj), and lattice data j2()| . 



The absolute value of the Fourier transform of the (transverse) gluon propa- 
gator with respect to Euclidean time is shown in Fig.|2| First, we clearly observe 
positivity violations in the gluon propagator. Second, the agreement of the nu- 
merical Schwinger function with the Fourier transformed fit is excellent. The 
crucial property of the gluon propagator is that it goes to zero for vanishing mo- 
mentum. This can be seen from the relation = D{p = 0) = / df^x D{x) (with 
D(ji) = Z(ji'^)/p'^) which implies that a nontrivial propagator function, D{x), in 
coordinate space must contain positive as well as negative norm contributions. 

The function Q contains only four parameters: the overall magnitude which 
due to renormalization properties is arbitrary (it is determined via the choice 
of the renormalization scale), the scale AqcD, the infrared exponent n and the 
anomalous dimension of the gluon 7. The latter two are not free parameters: k 
is determined from the infrared properties of the DSEs and for 7 its one-loop 
value is used. Thus we have found a parameterization of the gluon propagator 
which has effectively only one parameter, the scale AqcD- 



5 Spectral Properties of the Quark Propagator 

When discussing the results for the quark propagator we stated already that 
a dressed quark-gluon vertex is mandatory. For the present discussion it is 
important to note that such or similar solutions of the Slavnov-Taylor identity 
for the quark-gluon vertex will always result in the appearance of a quark-gluon 
coupling term A_B proportional to the sum of quark momenta, 

Vf^\p,q) :=SA, + AB,-f ... 

= ^ 7>^ + « 2 2 (p + 9)^ + .... (8) 

2 p'^ — 

Such a coupling, being effectively scalar, is per se not invariant under chiral 
transformations contrary to the leading term, EA, of the quark-gluon vertex. 
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Figure 3: The Fourier transform of the scalar part of the quark propagator for 
three different types of vertices. Positivity violations disappear once the scalar 
coupling, Ai3, is included in the vertex |21j . 

It is important to realize that the term Ai? appears only in case chiral symme- 
try is already dynamically broken. Thus it is consistent with the chiral Ward 
identities. Its existence, on the other hand, provides a significant amount of 
self-consistent enhancement of dynamical chiral symmetry breaking. Fairly in- 
dependently of the form of the gluon propagator the resulting quark propagator 
respects positivity if the term AB is included in the quark-gluon vertex |21|. 

A peculiar feature of the Fourier transform of the scalar part of the quark 
propagator is the curvature appearing for small Euclidean times. There are at 
least three possible sources for this. A leading singularity on the real momentum 
axis may be accompanied by additional real singularities at larger masses or by 
complex conjugate singularities with a larger real part of the mass, or it may 
be the starting point of a branch cut on the negative real momentum axis |21| . 

6 Extrapolation to the chiral limit 

Lattice simulations of dynamical chiral symmetry breaking are facing problems 
from two sides. Finite volume effects tend to obfuscate results in the very in- 
frared. Furthermore, although in principle the Ginsparg- Wilson relation enables 
one to implement chiral quarks on a lattice, small quark masses are computa- 
tionally very expensive. Therefore lattice simulations are usually carried out at 
finite bare quark masses and then linearly extrapolated to the chiral limit. The 
DSE-approach, however, suggests that this linear extrapolation is inaccurate 
\'22\. The DSE for the quark propagator also contains the nontrivial relation 
between the dynamically generated quark mass at small momenta and the renor- 
malized quark mass at a perturbative renormalization point. To illustrate this 
point we display the dynamical mass M(p^, to^) from the DSEs for two different 
momenta, — and — 0.38 GeV^ in Fig. 0| One clearly observes the cur- 
vature in the DSE-results for both momenta, reflecting the nonlinear behaviour 
of the underlying equations. We compare these results with lattice data taken 
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Figure 4: Relation between the dynamically generated quark mass at = 
0.38GeV^ and = to the renormalized quark mass at fi — 19 GeV from 
the DSEs and from the lattice data of ref. JH| 

from ref. [T^ and from ref. |^. The lattice data are consistent with a 
linear fit. However, by multiplying our DSE-results with an overall factor we 
obtain a curved mass function that mimics a DSE-model-fit to the lattice data. 
The extrapolated chiral dynamical quark mass at — 0.38 GeV^ is roughly 
20% below the one obtained from the linear fit. Employing only the leading 
7^-structure of the quark-giuon vertex, as done in ref. "22", leads to even more 
drastic deviations. Thus we confirmed that the linear extrapolation of lattice 
data to the chiral limit cannot be trusted well. 

7 Epilogue 

Recent years have seen a lot of progress in Strong QCD. In Landau gauge we 
gained an understanding of the infrared behaviour of gluon and quark propa- 
gators. We have verified the Kugo-Ojima scenario for gluon confinement and 
found dynamical chiral symmetry breaking in an ab initio calculation. 

We have proposed relatively simple functions to describe the running cou- 
pling, the gluon and the quark propagators (see ref. [2] for details) for all 
possible values of momenta. These have the potential to provide a basis for a 
hadron phenomenolgy based on quarks and gluons, even and especially in the 
non-perturbative regime. 
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